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ABSTRACT 

We exami ne the problem of modelling tunnelling phenomena in more than 
one dimension. We find that existing techniques are inadequate in a wide class 
of situations, due to their inability to deal with concurrent classical motion. 
We show how to generalise these methods to allow for complex momenta, and 
demonstrate the improved techniques with a selection of illustrative examples. 
We comment on possible applications. 
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Tunnelling is perhaps one of the more versatile concepts of quantum theory. It may 
be used in tunnelling microscopy 1 , or in describing a phase transition in the universe 2 ! 
But, whatever tunnelling problem one wishes to address, all have in common the idea 
that two regions of space classically separated by a potential barrier are not quantum 
mechanically isolated, provided the barrier separating the two is finite. The reason is that in 
the classically forbidden regime the probability density is non-zero, although exponentially 
damped. 

Such is the simple picture. In practice however, calculating a tunnelling amplitude 
accurately is not so easy, except in a range of special cases, essentially one-dimensional 
problems or problems which can be recast as such. Recently, in considering inflation- 
ary scenarios, people have been examining models containing not only an inflaton scalar 
field responsible for the false vacuum energy, but also an additional scalar field, either by 
extending the gravitational sector, as in extended 3 or hyperextended 4 inflation, or by in- 
cluding an extra scalar field, as in double field inflation 5 . These ideas have in common the 
notion of a single field whose tunnelling is influenced by the classical evolution of a second 
field. This allows for a time dependent nucleation rate which resolves several problems 
with the old inflationary models. To date, most calculations 4-7 have involved a “freezing 
out” of the second field, merely using it to provide dynamical evolution parameters in the 
one-dimensional problem. Since most tunnelling calculations involve stationary state Eu- 
clidean time techniques, these more complex multi-dimensional models beg the question 
as to whether such techniques are really valid when there is classical evolution. The aim 
of this paper is to approach this problem in a way that treats both types of evolution on 
an equal footing. 

In this paper, we explore this question for the test case of quantum mechanical tun- 
nelling, and consider a variety of illustrative potentials. We concentrate on the stationary 
phase or quasi-classical approximation, examining how, and under precisely what condi- 
tions this may be applied. In particular, we try to avoid assigning Euclidean time any 
preferred status, keeping it firmly in the category of an optional mathematical tool. In 
dealing with this more general class of problems, we find that we have to modify the WKB 
matching conditions to allow for the passage of real momentum under a barrier and com- 
plex momentum beyond. This naturally increases the complexity of the process of solving 
for the wave function under the barrier and beyond. We develop a method for solving 
these two problems and apply it to a variety of two-dimensional examples. 

The layout of the paper is as follows. We begin by reviewing the quasi-classical 
approximation, and what limitations it places on the type of tunnelling wave functions to 
be considered. We also review the role of Euclidean time, and highlight the problems of 
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trying to apply these techniques to more general multi- dimensional barriers. This leads us 
to our description of the identification of turning points in WKB solutions and derivation 
of appropriate matching conditions. We then solve the Schrodinger equation under the 
barrier in the stationary phase approximation. Finally, we apply our techniques to a few 
simple, but illustrative, examples. 

1. Tunnelling and the Quasi-Classical Approximation. 

Before reviewing the quasi-classical approximation, we will recap on what is usually 
meant by “tunnelling”. Tunnelling amplitudes for scattering problems are reasonably well 
defined as the ratio of the amplitudes of the emergent and incident wave functions. For 
‘stationary’ problems (a particle ‘tunnelling’ from one well or channel to another) we will 
consider a particle to be tunnelling out of its initial well if the probability of finding it in 
that well or channel goes predominantly as e -rt /\ T is the tunnelling amplitude. 

Now let us review the quasi-classical approximation®. As the name suggests, this ap- 
proximation extracts the leading order “classical” behaviour of the system, and is only valid 
when the de-Broglie wavelengths of the particles are small compared to the characteristic 
scales of the motion. Explicitly, if we set 

if} = e i"lh e -iEt/h ( 1 . 1 ) 


the Schrodinger equation gives 

( 1 . 2 ) 

as the equation of motion for <r. The quasi-classical approximation is that we drop the 0(h) 
term in (1.2), which clearly requires that | V<r| 2 » h\ W|. If we identify the momentum p 
with Vcr, we see automatically that the approximation is invalid for very small momenta, 
i.e. in regions where E ~ U , Thus a classical particle at rest has no stationary quasi- 
classical counterpart. However, what is perhaps less obvious is that this approximation 
also breaks down for more general turning points, when the momentum becomes parallel 
to the potential. We will discuss this point in more detail later. 

In general, if we are given <T{ at some initial point x,-, we may write the solution for a 
as 

<r(x)= f Vfr.dl + <t { (1.3) 

where dl is any path interpolating between X{ and x. This is all very well, but what is V<r? 
The Schrodinger equation is a scalar equation, relating the magnitude of Vcr to £/, yet we 
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presumably need to solve it throughout some multi-dimensional region of space. Banks, 
Bender and Wu 9 solved this problem for the case of tunnelling from a localised state in 
more than one dimension. We will summarise their approach (and others 10 ) here, before 
going on to discuss tunnelling from non-localised states. 

Essentially, they treated (1.2) rather like a geodesic problem, setting V<r oc dl, thus 
making (1.3) a scalar integral with respect to the length parameter along the path. This 
requires that Vcr has constant phase, and hence is purely real or imaginary, the latter case 
corresponding to tunnelling. To solve this tunnelling problem, they set er = ig and thus 
obtain to lowest order in h 


(Vg) 2 = 2m{XJ — E) 



^2m(U - E)dl, 


(1.4a) 
- (1.46) 


where 1 is the path which minimises the integral, the escape path , and x^ is the point 
of emergence of that path into the classical regime. In this case, it is clear how to solve 
(1.4a) exactly, however a common alternate method, indispensible in field theory, uses a 
Euclidean time description. 

First we replace p e = Vg in (1.4a), which yields 


H = 


2m 


p l + U = E. 


(1.5) 


'This can then be interpreted as a Hamiltonian problem of particle motion in the (inverted) 
potential —U. We then see that there is a non-trivial solution, x(r), which interpolates 
between the initial position of the particle and its position of emergence from the barrier. 
With our Hamiltonian interpretation we identify x(r) with p c , and r is thought of as 
a time parameter. The identification of r as a Euclidean time comes from noting that 
= p = ip e = tx(r), thus r = it. Solving the Hamiltonian problem along the trajectory 

gives: 

+ U = E (1.6)- 


hence 




p.xdr = i J 2(U — E)di 


(1.7) 


The advantage of this approach is that it provides a straightforward means of calculat- 
ing the escape path and action in terms of a classical mechanics problem, the disadvantage 
is that it introduces a fictitious time parameter along that path which can introduce con- 
ceptual confusion in trying to interpret what the ‘particle’ is ‘doing’. For example, a 


4 


statement often made is that because Euclidean time is somehow orthogonal to real time, 
tunnelling happens instantaneously; this is rather confusing. Tunnelling is a statement on 
the dynamical evolution of the wave function according to quantum mechanics, whereas to 
say “happens instantaneously” suggests an observation, i.e. an interaction with the system. 
The only process one can label as tunnelling is the quantum mechanical leakage of proba- 
bility across a barrier which is distinct from the path of a classical particle, and also distinct 
from an observation of the particle on the other side. The actual question of tunnelling 
time is an extremely subtle and complicated one (for a review see Hauge and Stovneng 11 ), 
and to some extent depends on how one chooses to formulate it. However, what one can 
say is that a probability density at time fo -b e a-t some point in space other than that at 
which a particle was observed at time t 0 is not a demonstration of instantaneous motion, 
but rather a reflection of the dynamics of the Schrodinger equation. 

Before going on to describe a more general approach to ‘escape paths we will make 
the obvious remark (which Goncharev and Linde* 2 discuss in more detail) that the stan- 
dard Banks- Bender- Wu method requires that (V^r) 2 = x 2 — > 0 at each end of the escape 
path, otherwise the solution cannot be matched with a classical ‘real time solution in the 
asymptotic regime. In the case of escape from a localised state, this can be guaranteed, 
since one end of the escape path is necessarily fixed and the other end is varied freely to 
find the minimum action, therefore we are free to set x 2 — 0 at each end. However, in a 
more general scenario, such as channel-channel tunnelling, this need not be the case. Both 
ends of the escape path can now be moved, and hence it requires two initial conditions 
(rather than one as in the previous case) to fix the start of the path. Thus, since our 
equations of motion are second order, we have used up our quota of boundary conditions, 
and we just have to hope that in varying our other endpoint we hit lucky . Unfortunately, 
from (1.6) we see that x — VC/ along the path, thus if U is for example monotonically 
increasing parallel to the barrier, x has a strictly positive component in that direction. 
Therefore x 2 can never be zero at both ends of any trajectory and the Banks- Bender- Wu 
method is not applicable. 

Another disadvantage of this approach is that it makes no allowance for the transport 
of real momentum under a barrier. In the case of a continuous symmetry it is clear that such 
transport does occur, since the extra degrees of freedom decouple from the problem, and 
the wave function is merely a product of a tunnelling wave function with suitable transverse 
momenta eigenfunctions. However, trying to incorporate this into a fully Euclidean picture 
is not so easy 13 . And of course we have totally ignored the problem of picking an initial 
position from which to integrate, since we can only localise a particle at the expense of 
information about its momentum. We clearly need a more general picture of tunnelling. 
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As we have noted above, in the approach of Banks, Bender and Wu, we must take space 
to be divided into regions in which Vcr is entirely real or entirely imaginary; this is what 
allows us to make the identification of V<7 with some x(“time”). Clearly we then need x 
to vanish on the boundaries where V<7 changes from real to imaginary, but this obviously 
means that we have restricted ourselves to a certain subset of problems. As we have seen, 
it is not difficult to find a potential which does not fit into this subset. We must find 
an approach which allows for complex momentum as well as dealing with the problem of 
matching between classical and tunnelling regimes. This is what we will now develop. 

2. Tunnelling with classical motion - general formalism. 

We now turn to how we can modify the escape path techniques for more general poten- 
tials. The class of potentials we will be interested in are ones which contain two channels 
(or asymptotic regions) separated by some barrier in the ‘x’-direction. We envisage that 
this barrier has some arbitrary y-dependence, but that this is secondary to the height of 
the barrier. Our main assumption, other than that of quasi-classicality, is that this barrier 
is always sufficiently high so that the division between ‘classical’ and ‘tunnelling’ motion 
is dear, i.e. U E. 

Recall that the essential problem of tunnelling with classical motion is twofold. First 
the transition from classically allowed to classically forbidden regimes is no longer char- 
acterised by real momentum becoming imaginary (or E ~ U): the momentum retains 
both red and imaginary pieces. Secondly, because the ‘momentum’ is now complex, the 
evolution equations become more complicated. The first point, deciding on matching con- 
ditions, is crucial; we may fix (Vy) 2 = 0 as we enter the barrier, but we must be able to 
interpret emergent solutions with non-zero Vy. 

In order to deal with complex momentum in the quasi-classical approximation, we 
rewrite (1.2) in terms of the real variables / and y, where a — f + ig: 

(V/) 2 - (Vy) 2 + fcV 2 y = 2 m(E - U) (2.1a) 

2V/.Vy — fi.V 2 / = 0, (2.15) 

this latter equation representing the constraint that the probability current, j = e~ 29 ^ V/ , 
is divergence free. 

We now turn to the problem of matching conditions: where they should be applied, and 
how to determine them. First let us describe how to identify a turning point. Physically, 
the breakdown of the quasi-classical condition occurs when the de-Broglie wavelength of 
a particle becomes comparable to the physical scales of the potential. In problems where 
the momentum changes from being totally real to totally imaginary (or vice versa), this 
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is characterised by E ~ U. However, if we wish to allow' simultaneous ‘classical’ evolution 
with tunnelling (i.e. complex momenta) then we must be more specific about the quantum 
to classical transition. 

In the ‘classical’ regime we expect the variation of the phase of the wave function 
to do mi nate, whereas in the tunnelling regime the variation of the amplitude should be 
dominant. In each of these regimes, the quasi-classical condition breaks down when 

|V<r| 2 ~ W| (2.2 a) 


in other words when 

(V/) 2 ~ hV 2 f 
(Vg) 2 ~ KV 2 g . 


( 2 . 26 ) 


These relations are satisfied either when |V/| or |V^| are small, i.e. when E ~ U, or when 
|V 2 /| or |V 2 flf| become large, this latter situation occurring near generic turning points 
of the motion, approximately when the momentum becomes essentially orthogonal to the 
gradient of the potential. To see why this is, let us consider an incoming plane wave 
scattering off some barrier. The integral curves of V/ trace out the path that a classical 
particle would follow in that potential. The breakdown of the approximation occurs when 
neighbouring, initially parallel, trajectories cross. (These can be either /-lines or < 7 -lines.) 
The set of points where rays cross each other is referred to as a caustic and it is at these 
points that the approximation breaks down due to V 2 / (or V 2 <j) becoming unbounded. 
(In the path integral language, as we pass through the caustic the path changes from a 
minimum of the action, 5, to a saddle point; that is one of the eigenvalues of 6 2 S becomes 
negative. In the expression for the propagator there is a term of the form (det 6 2 S) 3 
which diverges 14 .) The matching conditions should strictly be applied along the caustics 
of the / or g curves. However, since the barrier is assumed to be slowly varying along its 
length compared with its steepness, this actually coincides with the turning point of the 
/-line, when the momentum is orthogonal to VZ7. 

The matching conditions in the situation E ~ U are the standard matching conditions, 
that is, our initial conditions for integrating away from the turning point x<j are g = 0, 
V <7 = 0, V/ = 0 if matching from classical to tunnelling, and g = g(x/), Vg = 0, V/ = 0, 
if matching from tunnelling to classical at X/. In the second case, where there is some 
transverse momentum, then continuity of the wave function demands that the momentum 
parallel to the barrier is preserved. If this momentum is entirely real (or imaginary) the 
orthogonal part is matched as before. More precisely, the wave-function is given locally by 


i p = e Ai ( h s (2 mU, x )»z) , 


(2.3) 
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where Ai(a:) is the Airy function. This gives an extra phase of e 1 ? between the incident and 
reflected wave, which is a generic feature at such turning points and can be traced back to 
the change of phase in the determinant factor of the propagator caused by reversing one 
of the eigenvalues of 6 2 S, 

In the case where the incident wave has complex parallel momentum, a y 4- ib y , we 
might expect that 

v> = e^+’^Ai (h~ 2 ' 3 (2mU, x )*/*(* + + ( 2 . 4 ) 

However, using the asymptotic expansion for the Airy function, we find that the momentum 
of ip in the x-direction is given by the energy equation 

P 2 x + ( a y + iby) 2 = 2 m(E - J7(x)) , (2.5) 

hence if a yy b y ^ 0 p x can never vanish. Returning to the functions / and g in equation 
(2.1a), we see that whenever the distance between the /- curves scales as d along their 
length, Vg scales as d” 1 , thus as we approach a ‘caustic’ of /-lines, d —► 0 and Vg 
diverges, unless it happens to be zero. We can therefore think of the (V^) 2 term as 
a repulsive potential between the /-lines, which provides a smooth transition between 
different regimes of our approximation. 

In the examples that we consider the walls of the barrier are step-functions. In this 
simplified situation, we fix the (possibly complex) parallel momentum, and then compute 
the momentum perpendicular to the barrier using equations (2.1a,b), choosing the solu- 
tion corresponding to an exponentially decaying wave-function under the barrier and an 
outgoing wave at the far side. In both cases, if we label the incoming and outgoing perpen- 
dicular component of the (complex) momentum by k\ and respectively, then a trivial 
calculation matching the wave-function and Its derivative at the boundary shows that the 
reflected and transmitted waves have coefficients relative to the incident wave of 

ki — fco 2fcj 
Jci + Jfeo ’ h + k 0 

respectively. Normally, the wave reflected from the incident decaying wave under the bar- 
rier is only an exponentially small correction (associated with multi-instanton solutions) to 
the wave-function, and so is generally disregarded in what follows, although it is responsible 
for the very small but non-zero “tunnelling component” of the probability current. 

Now that we have initial conditions for Vg and V/ under the barrier, and also at the 
far edge of the barrier, let us turn to the equations of motion. 


( 2 . 6 ) 
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Since we are assuming E <C U, Vg is clearly dominant under the barrier, however, 
this does not mean that we can neglect V/ as an 0(h) correction, otherwise we would use 
existing techniques. Instead we want to consider a situation where | < 1, therefore 
we adopt a step by step procedure in solving 

(V<7) 2 =2mU -{2mE -(Vf) 2 ) (2.7 a) 

Vf.Vg = 0, (2.76) 

bearing in mind that (V/) 2 , E are of the same order, and small compared with (V^) 2 , U. 
The first step is to find the leading behaviour, that is, to solve 

(Vgf = 2 mU . (2.8) 

This is solved using the existing techniques: we use the momentum transfer equation 

V Vg Vg = mW (2.9) 


to find the integral curves of V<7, and then integrate the scalar g along them: 


9 = 



( 2 . 10 ) 


This gives us the solution to leading order. The next step is to use the initial conditions for 
V/ to integrate (2,7b) through the barrier, since (2.7b) merely tells us that / is constant 
along (j-lines: / = f(<p) where y 0 = ¥>(x,y) are the integral curves of g . Finally, we input 
this solution for / back in to (2.7a) to obtain the correct form of g to order E. 

Now that we have a systematic method for solving underneath the barrier, an ap- 
parently reasonable question to ask would be what the actual flows of the particles were 
across the barrier. This, as it turns out, is a very difficult question to address, since it 
tries to relate a classical notion (a path) to a classically forbidden regime. Even setting 
aside such interpretive reservations, we see that our solution as it stands cannot represent 
tunnelling, since the probability current, j = e 2p /*V/, is orthogonal to the integral curves 
of g and can therefore never leave the turning point. Clearly our solution is incomplete, 
since we know that particles do in fact tunnel, albeit with a very small amplitude. Let us 
first highlight a ‘trick 5 by which we may obtain the form of the probability current before 
discussing its true origin and therefore limited interpretational value. 

Consider first (2.1b) in one dimension. If we were to expand / in powers of h we 
would conclude that / = 0. However if we directly solve (2.1b), we obtain 


S't = &(*/)«’<''" ,/l • 


(2.n) 
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Thus, we find that as well as the zero solution, we also have a non-zero g-dependent phase 
under the barrier. This suggests that in general the probability current along a (y-line is 

j T oc Ke-***M h (2.12) 


a constant. We also have an order of magnitude estimate of the phase change along a 
g-line as hK, where K is the outgoing momentum. Note that all the phase change takes 
place at the boundaries of the barrier, therefore to get the exact phase change we should 
look at the matching conditions. In the case of a step function, (2.6) implies 


A/r = 0 = arg 


2kj 

ho + hi 


2ki 1 

N ho + h f F j 


(2.13) 


where the subscripts N and F refer to the near and far sides of the barrier respectively. An 
obvious point which we will nonetheless make is that if a path does not emerge from under 
a barrier, g increases without limit along that path, therefore the appropriate solution for 
f T is the zero solution. This implies that no tunnelling can occur along such a path. 

It would be nice to associate this probability current with the path of the particles, 
however, this would be hopelessly incorrect. We can only associate the probability current 
with a path if we have a single WKB-like wave function. Once we have a superposition 
of wave functions, such as incoming and outgoing waves, the probability current only 
represents a nett flow due to the interference between the various waves. In the case 
of real momenta, we can directly sum the probability current associated with each wave 
function to get the nett current. With complex momenta however, which unfortunately 
is precisely what we are interested in, it is the interference which gives the probability 
current. Thus, although this trick which we have illustrated does extract a probability 
current from a complex momentum wave function, the non- vanishing of this current relies 
crucially on the existence of a point of emergence for this wave function into the classical 
regime, and this is exactly equivalent to the existence of another, exponentially growing, 
branch of the solution. As such V/r cannot be interpreted as a flow of particles under the 
barrier, but rather as a shorthand way of estimating the multi-instanton corrections which 
does in fact give the correct order of magnitude for the tunnelling current (although it is 
out by a factor of order unity). 

Finally, before investigating a few concrete examples, we should summarise the limits 
to our approximation. First, we have assumed E <C U in order to facilitate the solution 
of (2.1) as an expansion in E/U. Our second main assumption involves the matching 
conditions in terms of the local orthogonal coordinates along the caustic, this requires that 
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the caustic not be too strongly curved which translates into a restriction on U yyy . From 
(2.4) we see that the quasi-classical approximation is invalid for |x| < h 2 ^(2mU, x ) -1 / 3 . 
This gives a rough order of magnitude limit U, yy < ft~ 4 / 3 (2mJ7, x ) 2 / 3 , which one can sub- 
stantiate with a more careful calculation. This shows that our turning point treatment is 
generically valid. In practise, calculational complexity will be the limiting factor for Z7, vy 
- curved caustics are more difficult to deal with. 

3. Plane wave scattering - examples. 

In order to examine the flows of the wave function we first consider a steady state flux 
of particles impinging on a variety of barriers, that is, We do this because 

we no longer have localised quasi-classical states in the direction of classical motion. We 
follow the wave function under the barrier and out into the second ‘classical’ regime. After 
calculating the form of the tunnelled wave function, we use the momentum eigenstates to 
build the physically more realistic situation of a Gaussian wave packet hitting the barrier. 
This more realistic scenario allows us to highlight the existing controversy of tunnelling 
times for the square barrier, as well as illustrating some interesting new properties of more 
general barriers. 

To illustrate our techniques we have chosen three examples in order of increasing 
complexity. We begin with the simplest possible two-dimensional potential - a separable 
square barrier. This allows us to check our calculations against an exact solution, and 
also leads naturally to the second example: a square barrier of varying width. Finally, we 
consider a square barrier of fixed width and varying height. These latter two potentials 
demonstrate some very peculiar scattering properties as we shall see. 

Example 1: J7(x) = V r 0(x)©(a — x). 

For this potential we know the exact form of the stationary eigenfunctions: 

i>(x) = e ip * v/h xl> Pi (x) (3.1) 

where 

{ e ip x z/h _ e 'p x a l h e~ x *\ x f h sinh/ca/h x < 0 

\Te ip ^ h [(1 + + (1 - *>/*] 0 < x < a t 3 ' 2 ) 

Te'Pi x / h x > a. 

is the usual one-dimensional square barrier wave function, with 
p x 2 = 2m E — Pj 2 k. 2 = 2mV — p x 2 

T — (cosh ac a/h -(— )sinhAca/ft 

V 2 K Pi 
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It will be convenient to rewrite this transmission coefficient as 


_ -ip 1 a/h -iQ/h 

D 


where 

D 2 = 4p, 2 /c 2 + 4 m 2 V 2 sinh 2 /co/ft 
gives some measure of the transition amplitude, and 


0 

h 


tan 1 



tanh 


/col 


(3.4) 


(3-5) 


(3.6) 


gives the 'transition phase’. This latter quantity is important in discussing some interpre- 
tations of the tunnelling time. 

For V > E we see that 


D ~ 2 mV sinh — ~ mV e Ka ^ h 

h 


0 

h 




tan 1 



* 2 Pi 
2 /c 


(3.7a) 

(3.76) 


and hence 

Comparing the transmitted solution with the incident solution, we see that the transmission 
amplitude, |T|, is given approximately by ~~^ e ~ Ka (having set k ~ s/2 mV) and that the 

transmitted wave acquires a phase ( " ■ % — — f + -£*-)• 

To solve the Schrodinger equation using a quasi-classical approach is in this case very 
straightforward since U is a function of one variable only. The solution of (2.7a) to leading 
order is g = \/2 mVx. Equation (2.7b) then implies that / = p,y. Equation (2.7a) then 
implies that g - */2m(V~^ E) + p, 2 ® as required. 

To solve in the asymptotic regime we note that since V<7 = 0 along x = a, g must be 
a constant, <7/, and therefore V/ returns to its original value (p,,p, ). To find the phase 
shift across the barrier, O, we use (2.13) and obtain A/r = 0 = tan -1 j in agreement 
with (3.7b). Thus, integrating out from x = a, we obtain 

g = Ka (3.9a) 

/ = p,y + Pi(x - a) + fi.0 (3.96) 


in agreement with the exact solution. 
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Now we will attempt to solve for potentials which have a small y-dependence included: 
Example 2: ^(x) = 70(i)0(a — x — ey). 

This represents a barrier of constant height and varying thickness. Clearly, since our 
initial conditions for integrating under the barrier are the same as for the first example, 
the solution for / and g under the barrier is the same: 


/ = P a 3/ i 


g = yj2m{U — E) + p* x = kx 


(3.10) 


To find the solution in the asymptotic region, we must recall the boundary conditions to 
be applied at the fax edge of the barrier: V/ |( , Vy (| are preserved. Now, if we assume that 
e 1, and only keep terms of order e, then the tangent and normal vectors to the fax 
surface axe, respectively, 

t =(T) n =(1)- (3U) 

Therefore 

Vy M = -cacT and V/ |( = p,T . (3.12) 

We may then use equations (2.7a,b) in the asymptotic region to conclude that 

/ou. =P 1 {x-a + ey)+ p,(y - e(x - a)) + he 

= Pi£ + P» T ? + ft© (3.13a) 

9o*< = - cut] + — e/c£ (3.136) 

Pi 


where 

£ = x — a + ey 

7j = y - e(x - a) 

and 0 = tan -1 2p ~* €p is given by (2.13). The solution for g in the asymptotic region 

is now a function of 77 and £ (or x and y) This is in contrast to the previous example of 
the square barrier. It is intuitively obvious that the the probability density decreases in 
those regions where the barrier is thicker. The spatial variation of g will lead to interesting 
phenomena when we come to discuss the scattering of a Gaussian wave-packet off this 
potential in the next section. 

Example 3: £^(x) = Ui(x)e~ ey . 

In this example we have a barrier of constant thickness, but varying height. The 
matching conditions for a step- function give the initial conditions for integrating under 


(3.14) 
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the barrier: 


V/| v=0 + = (0,p 2 ) 

Vg\ y=o+ =(V2mVe-^,0) 


(3.15) 




^( 0+ »J/) = e * 


2Pi 


lp t + i\/2mV e~^ 


The first step in solving equation (2.7a) is to find the integral curves of Vg. If we 
write x = Vp (cf Euclidean method) then by examining x = VI7 we see that 


y 2 = 2mV(e~ ty - e' y °) 
x 2 = 2 mVe~ fy ° . 


(3.16) 


The equation for the integral curves is therefore 

^ = -y/ e -<v-vo) _ i . (3.17) 

(3.18) 


this latter approximation being valid when « < 1. Note that the form of the integral 
curves is invariant along the barrier. The curves all asymptote the line ex = 7r, so we 
will impose a < 7 r/e to ensure that (V/) 2 remains small. Now, for these curves, = 
e -*(y 0 -y)/2 cos *± m Hence 

g= ! VWds = / s/2mVe~ ty / 2 sec y <£x 

J Xq J Xq ^ 

= \j2mV e “ c ^°/ 2 — tan — . 

e 2 

= V2rnV e~ €y f 2 — sin — . (3.19) 

e 2 

Having found the zeroth order solution for we must now apply (2.7b) to find /, 
using / = p 2 y Initially. From (2.7b) we see that / must be constant along integral curves 
of 5 , hence 

2 ex 

/ = pAy - - lo s cos y) • (3.20). 


This has solution 


2 , ex 

y = yo + - log cos — 
e « 


- yo - 


ex* 
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From this we may deduce that (V/) 2 = p 2 sec 2 ~ and hence 


(V^) 2 = 2 mV e~ ty + p 2 sec 2 — — 2mE . 

4 


(3.21) 


Now we find ourselves in the situation of having a modified equation for the integral curves 
of Vg. 


dy 

dx 


(e-^-vo) - 1) 


\ 1 - f e*v° + j^e'vo sec 2 f 
1 + — e £yo 

* 1 Tf 


= — ^/ e -<(v-Vo) _ 1 




, e'* u sec" ^ 
4mV 2 


Setting 




we transform (3.22) back into (3.17), with x replaced by X . Hence 


2, tX 

y = yo + - log cos — 

C 4 


are the new integral curves, and 

g = f \Pl TnVe~ €y °/ 2 sec 

Jx o * Y 


p 2 

1 + -^-e f vo 

2m7 


= £ «•* f 0 + Sv'™) dx 

= v^e-/= sin f (l + ^e*» sec 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


is the new <7. 

Having now verified the stability of our solution, for calculations! brevity we will now 
retain only the leading order parts of / and g in what follows. We must now match the 
wave-function across x = a. We know that /, g, all match continuously across the 

boundary. In particular 


dy 

*2. = -V2^Ve~?sin- 
dy 2 


( 3 . 26 ) 
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at x = a + . Using our matching conditions, if we let p x = (|£ + iff) then 


p 2 = 2i? — (p, — iVZmV e ^ sin—) 2 

2 

= p 2 + 2mUe- ey sin 2 ^ + 2zp, V^Usin^e”? 

z z 


(3.27) 


There are two limits in which information can readily be extracted: 

( a ) yf§ > ~ 

In this case the barrier is very slowly varying with respect to the other scales in the 
problem (although not necessarily with respect to h), and the asymptotic regime is clearly 
identified as (Vp) 2 -C (V/) 2 . Clearly we may take sin^f ~ so that to order we 

have 

P, = Pi + *^v/2mTe-? (3.28) 

2 Pi 

Solving for / and g in the asymptotic region can be achieved by following the same 
procedure as under the barrier. In this case we first solve for f and then demand that g is 
constant along the integral curves of V /. This gives that 

2 to* 

f = p t (x-a) + p 1 (y- - log cosy) (3.29a) 

g = v'2mVae"^ [s_ ^ (l_a)1 . (3.296) 

These expressions are valid in the region 1 e[y — j^-(® — a)] log e. 

(b) ^ « Anf 

In this case we have 

p„ « V2 mV e~ ¥ sin^ + zp, (3.30) 

z 

In this limit one can ignore the energy term in comparison with (V/) 2 ,(Vp) 2 and the 
WKB equations for /, g become equivalent to the Cauchy- Riemann equations. We can use 
this to find an approximate solution for / and g in the region x > a. If g = g(a , y) and 

/ = /(a, y) on the boundary then letting z — y + i(x — a) we may write 


2 j . ta tlx — a) , 

g = Re^(a, z) + Im/(a, z) = -v 27711^11— cos e 2 + p 2 [x — a) 

6 Z Z 

2 / ta tlx — a) . 2 ta N 

/ = -Imj(a,z) + Re/(a, z) = - V2mVsiii— sin e 2 + p 3 (y log cos — ) 

t 2 2 t 2 


(3.31) 
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At this point let us make a few remarks on the form of this solution. Near the 
barrier the wave-function dies off doubly exponentially rapidly as y — > — oo, where the 
barrier is highest. This is not surprising since the barrier is exponentially rapidly growing 
here. However what is interesting is that the solution does not seem to make sense for 
x — » a + 7r / e as g — ► 0 here, making the wave-function (doubly) exponentially enhanced 
relative to its value near the barrier. This can be understood by looking at the integral 
curves of V/, which asymptote x = a + 7r/e. As we approach this line, we cross integral 
curves of V/ which emanated from x = a at larger and larger values of y, in particular 
where y/E/U ~ ea, where approximation (b) is no longer valid. Thus the solution (3.31) 
cannot be extrapolated to x ~ a + 7r/e, and can only be regarded as a solution close to 
the barrier. However, it is interesting in that it deals with a regime in which the barrier is 
rapidly varying in height. 


4. Scattering of Gaussian wave-packets 

We are interested in understanding what happens when a classical particle scatters 
off the sort of potentials we have been considering. So far, we have given a solution for 
incoming plane waves scattering off such potentials. The linearity of Schrodinger’s equation 
enables us to add these solutions together in order to find time-dependent solutions. We 
will be interested in considering solutions for which the incoming wave is a Gaussian 
wave-packet, since this can be thought of as representing a classical particle. After some 
introductory remarks we shall examine the solution for a Gaussian wave-packet scattering 
off each of the potentials considered in the last section. 

Briefly, we review Gaussian wave-packets to establish notation. Let us consider a 
solution to the one-dimensional Schrodinger equation, which at time i = 0 has the form 

'f(a:,0) = A«7 t“* e * e * (4-1) 


To find ^(ijt) we first take the Fourier transform 

F(k) = — [ dx 
2tt J 


= 2“’7t“Sa - * 


exp 




Then we have that 


ifcg till 

)e * e 


*(*,<)= j 

{ (x — 

- L — t: 
4n ^ 


(4.2) 


) 2 ip\t ixpi 


(4.3) 


2 m ti 


+ 
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where A = 2 ^ } . This solution represents a Gaussian w'ave-packet with momentum 

pi. The wave-packet has < Ax 2 >= 2ft 4 A|.A| 2 and < A k 2 >= ft 2 . The modulus of A is 
smallest at t = 0 and at this time the wave-packet is minimal. The wave-packet begins to 
spread appreciably for |i| > 

People are generally interested in two types of wave packet: those extremely peaked in 
momentum space ( A of order unity) and those equally spread in position and momentum 
space (A = 0(ft -1 )). Notice that even in this latter category, the wave packets are still 
very sharply (~ ft) peaked in momentum space. In the first case, the packet, although 
comparatively diffuse in x-space, maintains its shape for t < ft -1 , whereas in the second 
case, the packet starts to spread for t of order unity (although the spread is still of order 
ft); this is quite sufficient. We will therefore assume that Aft is at most of order unity. 

In the previous section we obtained the solution for a plane wave scattering qff three 
different potentials, what is now required is that we take the appropriate linear combina- 
tions of these solutions so that the composite describes an incident wave-packet which can 
be interpreted as a classical particle. This is simply a product of the solutions described 
above; that is 


*in (*,!/, t) = 


2h 2 A\fir\ 




exp 


{ 


4ft 2 A 


m m 


+ 


txpi iyp 2 

- l 


} 


(4.4) 

where E = ^(p 2 +p\)- This describes a wave packet with momentum (pi,P 2 ) hitting the 
origin at < = 0. The momentum profile of the wave-function has the form 


F(hi hi) - — l— — 

ft ft 27rft 2 \/jrA 


(4.5) 


In order to find the emergent wave-function we integrate 


* E (x,y,t) = JJ dfc 1 dfc 2 F(^- l ^)^out(fc il *,y> 




(4.6) 


where 


^out oc e* (/ou ‘ + ‘ 3oa,) (4.7) 

is the approximate form of the outgoing wave function calculated in section three. We may 
write this integral in the form 




27T h 


hrxJJ dhik 


(4.8) 
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where 

A(®, y,k i , fcj , <) = i/out Pout — 2\fi P* ) ~2rn^ i * (^*®) 

We use a saddle point method to calculate this integral making the exponent stationary 
with respect to Jfc t , k t . We are interested in the trajectory of the peak of the wave packet, 
the position of which is given by VReA = 0. We then use the set of equations 



VReA = 0 


(4.10a) 

(4.106) 


to determine the trajectory of the outgoing wave packet. 

Having now outlined the general procedure, let us return to our three examples, in- 
vestigating the behaviour of a Gaussian wave packet as it hits each in turn. 

Example 1: U(x) = V0(i)0(a — x). 

Here we substitute in (4.6) the outgoing wave function from equation (3.2): 


1/w = 


(4.11) 


where T is the transmission coefficient in eq(3.3b). Integrating over k 3 merely inverts the 
Fourier transform in the y-clirection since T is independent of y and the integral over the 
z-momentum is peaked about some k x close to p\ «C U. We therefore approximate T by 




2fc,<c(fc.) f ik t a ic(fti) i0(fci)\ , A 

„r - ^1"* ft ft - }' (4 - 12) 


This now simplifies the expression (4.9) for the exponent: 

A (l.l.ftj) =ik,(x - a) - - pi) 2 - Atf -rea-i0(fti) 


(4.13) 


Using the saddle-point approximation we expect that 


^ E {x,y,t) oc exp 


{■ 


Ah 2 A 


(y _^) 2 


+ 


*P iy 


771 


ip 2 t ) 

— } e n 
2mh J 


(4.14) 


where 


0 - 


dA 

dk x 




= i( x - a) - — (&* - pi) - 


ik x t 


+ 


ak. 


171 yJlmV — k x 2 


- iG'(k x ) (4.15) 
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In order to find the trajectory of the peak of the emergent wave packet, we combine this 
information with (4.10b) (which simply reduces to Imfc^ = 0). Doing this, we find that 
along the peak 

s = a + ©' + ^ (4.16) 

m 

where k ' x , the momentum of the peak, is determined by the real part of (4.11), which we 
may expand to order y/EjV to obtain 

k=p '{ 1 + tB ?) (4 ’ 17) 

This, together with the equation giving the peak in the y-direction which is obviously 
y = suggests that the peak emerges at time t e = —Q’m/k^ =0 (ft) at (a, and 

travels subsequently in a straight line. 

We have been deliberately obscure about the actual value of 0' (and hence t e ) since 
it would be misleading to advertise this as the actual tunnelling time when there is consid- 
erable debate on this topic 11 , indeed, as we shall see in the next two examples, taking this 
quantity seriously as a tunnelling time would lead to some interesting physical dilemmas. 
This time is called the phase time 15,18 of the tunnelling process, and is really only well 
defined asymptotically, well after the scattering/tunnelling process has been completed. 
One way of seeing why this must be so is to recall that the distortions of the wave packet 
close to the barrier make ‘definitions’ such as ‘when the particle hits the barrier’ subjective 
at best, meaningless at worst. Our reason for including this term is twofold. First we wish 
to emphasise the physical, quantum mechanical nature of tunnelling. Secondly, in the next 
two examples, we find factors in t e due to the variation of the potential and we would like 
to compare the two factors. 

Example 2: ^(x) = V0(i)0(a — x + ey). 

Having worked through the first example in some detail, we now summarise the steps 
for this example. Using (3.13), we may write the exponent (4.9) as 

A(£,J?,fc t , *:,,<) = - Ka + cKTj - i2-eic{ + »(fc t € + M + ©A). (4.18) 

where £ and q are given by (3.15). Making this stationary with respect to the yields 
the two (complex) equations: 

+ i(£ + 0,i h) = 0 

(4.19) 

-yr-(fc y -P,) - — K - ^ +i(v + ©,2 h) = 0 

A *h TTb fcj 


it 


-T r(*.-Px) 


■ k x T] 


kyK^ ky^ 


A h 


m 


K 


K 
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For the trajectory of the wave packet we use (4.10b) and (4.18) to find 


. K, K 

ImJb; = -eRe-£- 

K 

(4.20a) 

ImA:' — eRere . 

(4.206) 


Now we have a situation in which the momentum dominating the integral is complex. 
Writing k' z for ReJk^, and substituting from (4.20) for ImA:^, and similarly for k' y , the 
equations for the emergent wave peak can be seen to be 


4 = 


V = 


k' x t k'y „ ky 

»- t + t • 

(4.21a) 

h ± +i 

m 

(4.216) 


where 


e« 
A 1 


kL 


-Pi 



k' 

K y 


= P> 


(4.22a) 

(4.226) 

(4.22c) 


Thus we see that the packet emerges from under the barrier at a time and displacement 
given by 


m 

te = kL 


k 1 

0,i h 4- e-^(l — o.\K/k) 


k’ ( k ' 2 y 

Ve = Mir 0 ’ 1 -0 ’ 2 ) + M 1 + pH 1 ~ aA V«) 


subsequently travelling along the straight line 


LI t ( u '2 s 

V = ^ + i{l + ^(l-aXh/K) + 0(h) 


(4.23a) 

(4.236) 


(4.23c) 


with a damping relative to the incident Gaussian of 

e -*[o-«,.]/a ^ e -<tx ollt /A 


(4.23<f) 
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Tii ere are several interesting differences with the previous example. The first is that 
the t unnellin g time can now be negative! This peculiarity arises for a simple physical 
reason, the fact that the most energetically favourable time for the Gaussian to tunnel is not 
necessarily when the peak hits the barrier. Tunnelling amplitudes depend exponentially 
on the size of the barrier, therefore it is more favourable to tunnel when the barrier is 
thinner or lower. On the other hand, the probability density along x — 0 is damped by 
an exponential factor depending upon how far away that point is from the peak of the 
Gaussian. Clearly there will be a pay off between these two factors which may mean that 
it is more energetically favourable for the fringe of the Gaussian to tunnel, rather than its 
peak. 

The second thing to notice is that the wave packet emerges on the other side of 
the barrier at 77 = rj e , somewhat ‘downstream’ of where one might expect it. The two 
possibilities for emergence would be the perpendiculars from either the start of the barrier, 
or its end, that is, 77 = 0 or 77 = ae. Neither of these naive choices are equal to rj e . This 
shift has the same physical origin as the unusual t e . 

Finally, we should remark that the incoming and outgoing momentum are not quite 
parallel, since the (x,y) and (£,77) coordinate systems do not quite coincide. 

Example S: U^(x) = Ui(x)e~‘ v . 

We shall now discuss the problem of scattering a Gaussian wave-packet off a barrier 
of variable height. This problem is much closer to the sort found in inflationary models. 
As before, the emergent wave-function will be the superposition 

= JJ (4.24) 


where we include the factor of e so that the wave packet hits the barrier at (0,yo) 

rather than the origin. Recall ut was calculated in two separate limits. We now calculate 
the outgoing peak trajectory in each case. 

(a) » ea 

In this case, from (3.29), we find that A is equal to 


A = 


(k ( — P<) 2 


2Xh 


, i(k l (x-a)+k i (y—yo — - log cos ^))-V 2mVoe yo 

2 771 € £ 

(4.25) 
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Making this stationary with respect to the k { gives 


--^(fc» ~ Pl ^ ~ m k * + ^ X ~ ^ + T V2mV &( X “ a ) = 0 
- P.) -£*, + •(»- f 1°8 cos j) + jV2^V (l - ^) = 0 


As before, we use (4.10b) to find the trajectory of the peak of the wave packet 

Imfc; = ~^V2mVRe^ 

Z IC 9 

Imfc; = jV2mV 


Again, the momentum dominating the integral is complex. Writing fc; for 
before, the equations for the emergent wave peak can be seen to be 



2 , ea fc«t 

y = J/o + - log cos — + + e 

e 2 m 


where 


ea 


e = 


V2m7 

2XK 


K=?i 


eaXh 

T 


k' y =p 3 + ^V2mV 


Therefore the packet emerges from under the barrier at 

mefc' 


U = 


k'} 


2 ea / fc; 2 \ 

ye = yo + - log cos — + e I 1 + — 1 


subsequently travelling along the straight line 


2 - ea k' . » / fcy 

y = yo + - log cos — + — (x - a) + e I 1 + — 


(4.26) 

(4.27a) 
(4.276) 
4 etc. as 

(4.28a) 

(4.286) 

(4.29a) 

(4.296) 

(4.29c) 

(4.30a) 

(4.306) 

(4.30c) 
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with a damping relative to the incident Gaussian of 


e 


(4.30d) 


Thus, as in the previous example, the tunnelling time once more can be negative (for the 
same reason) and the peak of the transmitted wave packet emerges somewhat downstream 
from where we might have expected. 

In this case, from (3.31), we obtain 


2 — , £d .• «(? +<(■-•)) 2, ta ./_ _w (^» P») tk{i 

A = — V2 mV sin — e* ’ + tfc,(y -y 0 - - logcos — + t(x - o)) ^ — 

(4.31) 

Using saddle-point methods to evaluate this integral we arrive at the solution 


2, ta , _jx rz— —r . ta ( 1 e(x - o) t t . e(x - a)\ 

y = yo + - log cos — + e t y/2mV sin y ( — cos + — sin J 


x=a+ ^ +e ' !f ' / ^ Fsia f(-^ si 


e(x — a) t e(x — a) 

sin -i— - — - H cos — - 

2 m 2 


) 


(4.32) 


so that the peak emerges at x = a at a time and position given by 


te — 


XhV2mV sinf 


2 , ea ta 

y. = y. + 7 logcos T + — sm- 


(4.33a) 

(4.336) 


As in the previous examples there are three terms which contribute to the value of y at 
which the particle emerges: its initial position, the shift to this due to the curvature of the 
V</ integral curves under the barrier and a term which depends on the spread A. Again, 
the time at which the peak emerges can be positive or negative depending on the sign of 

P '*• 


5. Conclusions. 

In this paper we have generalised existing methods for calculating tunnelling processes 
to allow for complex momentum. This allows us to calculate tunnelling amplitudes in a 
wider class of potentials. We started by reviewing the BBW approach, then described 
how to generalise this to include complex momentum. The problem of solving under the 
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barrier was made tractable by assuming that the energy of the wave function was small 
compared with the height of the barrier. The matching conditions we derived by imposing 
continuity of the wave function allowed us to transport real momentum under the barrier 
and imaginary momentum beyond. 

In section 3 we applied this method to three examples: a square barrier for which 
we knew the exact solution, a step function of varying width, and then one of varying 
height. We found that the outgoing momentum was not necessarily parallel to the incoming 
momentum, a feature that we expect for generic barriers. Examining the scattering of 
Gaussians against these barriers allowed us to follow a wave packet tunnelling; probably a 
more realistic physical scenario. In recapping the square barrier case, we could illustrate 
the origin of one definition of tunnelling time, the phase time. We used this to explore the 
effect of barrier variation. This, as it turned out, was quite significant. Both in the case of 
the step barrier with varying height and that of varying width we found large corrections 
to t e and a shift in the place of emergence of the wave packet into the asymptotic regime. 
These turned out to have a simple physical origin. It is incorrect to assume that it is 
the peak of the incoming packet that dominates the emergent wave-function, and that 
tunnelling takes place when this peak is next to the barrier. Tunnelling is an exponentially 
suppressed phenomenon, and therefore it is far more favourable to tunnel where the barrier 
is smaller, even if the impinging probability density there is not a maximum. 

The last two examples we considered suggest that the definition of ‘phase time’ as 
it stands is not a good definition of tunnelling time since it assumes that the tunnelling 
process starts when the peak hits the barrier. It is possible that some modification of this 
description would give sensible results, compatible with the uncertainty principle for energy 
and time, although such a modification would naturally weaken the status of phase time. 
Unfortunately, we are not able to shed any light on whether some alternative definitions 
of tunnelling times 17-19 are any better. 

We restricted ourselves to the examples discussed since they illustrate the salient fea- 
tures of multi- dimensional tunnelling without too great calculational complexity. However, 
it was the cosmological applications of complex momentum tunnelling that originally inter- 
ested us, and these involve tunnelling processes from semi-localised states in field theory. 
The problem of calculating a tunnelling rate from a state localised in the x-direction should 
represent one step up in complexity from our examples and is currently being calculated. 
The problem of applying these ideas to field theory would involve translating the tech- 
niques presented here into the functional Schrodinger picture. This may prove to be very 
problematic, although clearly that is the next step in solving the two field problem. 
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